Radiative parton processes in perturbative QCD - an improved version of the Gunion 
and Bertsch cross section from comparisons to the exact result 
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In this work we compare the Gunion-Bertsch approximation of the leading order perturbative 
QCD radiation matrix element to the exact result. To this end, we revisit the derivation of the 
Gunion-Bertsch approximation as well as perform extensive numerical comparisons of the Gunion- 
Bertsch and the exact result. We find that when employing the matrix elements to obtain rates or 
cross sections from phase space integration, the amplitude by Gunion and Bertsch deviates from 
the correct result in characteristic regions of the phase space. We propose an improved version of 
the Gunion-Bertsch matrix element which agrees very well with the exact result in all phase space 
regions. 



I. INTRODUCTION 

When studying the dynamics of the partonic stage 
of relativistic heavy ion collisions based on microscopic 
perturbative QCD (pQCD) processes or on rate equa- 
tions for such processes, it has been found that elastic 
interactions alone are not sufficient to describe the evo- 
lution of the system Due to the numerical com- 
plexity of such dynamic approaches the incorporation of 
inelastic particle production and annihilation processes 
has mostly been based on leading order matrix elements. 
A commonly used approximation to the leading order 
pQCD matrix clement for partonic 2 f> 3 processes is a 
result derived by Gunion and Bertsch (GB) in 1981 
Q. This approximation gives a comparatively simple 
expression for the gluon radiation amplitude in terms 
of the transverse momentum of the radiated gluon k± 
and the transverse exchanged momentum q±. It has 
been widely used when solving transport problems (fj|-[l3l| , 
for example via rate equations or via microscopic trans- 
port approaches, most notably the microscopic transport 
model Boltzmann Approach to MultiParton Scatterings 
(BAMPS) dUni. In Ref. [HQ]! some efforts are made 
to go beyond the soft GB approximation while still ob- 
taining a relatively compact form. 

Recently, the results obtained within the BAMPS 
framework, especially with respect to the computed shear 
viscosity [2(| , have been challenged in a paper by Chen et 
al. [!l[ , who claim that a mis- or double-counting of sym- 
metry factors when applying the GB matrix element to 
inelastic processes might lead to an overestimation of in- 
teraction rates in BAMPS by a factor of 6. Also a recent 
work by Zhang [22j has addressed this issue and finds 
numerical discrepancies between cross sections based on 
the GB approximation compared to full leading order re- 
sults for the 2 — > 3 amplitudes, which have also been 
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known since the early 1980s [23|, |24[. The numerical dis- 
crepancy found in Ref. 1 221 is much less pronounced than 
the one claimed in Ref. [2l[. However, it is not fully clear 
how the results of these two works compare since both 
use different screening schemes. 

In this work we will address these issues in great detail 
by providing extensive numerical comparisons between 
the GB approximation and the exact leading order pQCD 
result as well as by analytically revisiting the derivation 
of the GB approximation. We will show that the approx- 
imation by Gunion and Bertsch needs to be carefully 
corrected in certain phase space regions when being em- 
ployed to the computation of cross sections or rates from 
the matrix element. To this end we propose an improved 
version of the GB matrix element which is valid in all 
regions of phase space. However, we want to empha- 
size that the deviations of rates computed in the original 
GB approximation compared to rates computed from the 
improved GB and the exact matrix element are caused 
deep within the approximations made by Gunion and 
Bertsch and are not given by simple symmetry factors 
as argued in Ref. (2lj . 

The paper is organized as follows. After a short revisit 
of the GB calculation we propose an improved version of 
the matrix element in Scc.|TTJ In Scc. lIIII the exact matrix 
elements for 2 — > 3 processes are introduced, whereas the 
phase space integration needed to obtain a total cross 
section is outlined in Sec. IIV1 Finally, we will compare 
the original and improved GB cross section with the exact 
result and discuss deviations. 



II. GUNION-BERTSCH MATRIX ELEMENT 

In this section we will argue that the approximation to 
the leading order pQCD matrix element for qq' — > qq' g 
processes (and more generally for other gluon radiation 
processes as well) as computed by Gunion and Bertsch 
[5[ needs to be carefully corrected when being applied to 
the computation of cross sections or rates. This correc- 
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Figure 1. Feynman diagrams for q + q — ► q + q' + g. 



tion consists of two parts, a) keeping a kinematic factor 
(1 — x) 2 , where x is the fraction of light cone momen- 
tum carried by the radiated gluon; and b) respecting the 
symmetry of the process by explicitly combining results 
from A + = and A~ = calculations, restricting the 
emission of gluons to the respective forward direction. 

We start by reviewing the derivation of the Gunion- 
Bertsch approximation, taking the example of qq' —> qq'g 
processes as in the original publication [f|. In leading 
order this process is given by the five Feynman diagrams 
shown in |Figurc l) 

In the following we label the 4-momcntum of the in- 
coming (outgoing) quarks with p\ and P2 (p3 and p±). 
The radiated gluon is denoted with k, while q is the in- 
ternal momentum transfer (the momentum of the inter- 
nal gluon propagator). Furthermore we are working in 
light cone coordinates where in the center of momentum 
frame the momenta of the incoming particles are given 
by 



Pi 



(Vi, 0,0,0) 

p%= (0,Vi,0,0) 



The Lorentz-invariant quantity 

x = k + Jp\ 



(1) 

(2) 



(3) 



characterizes the fraction of light cone momentum carried 
away by the radiated gluon. It can be related to the 
rapidity of the emitted gluon via 



k± 

7i e 



(4) 



The Gunion and Bertsch matrix element is derived 
in the high energy limit, i.e. the radiated gluon and the 



momentum transfer of the process are soft, 
kx < Vs 



(5) 



The third approximation explicitly stated by GB relates 
kx, qx, and x, 



xq± <C k± 



(6) 



With these approximations and from the constraint that 
all external particles are on-shell one can explicitly spec- 
ify k and q in light cone coordinates 



k» = ( x^Ts, P± - 



Xy/S 
^2 7,2 



k\jx + q\ - 2q ± kj 
(l-x)V^ 



(7) 
(8) 



The original computation of the GB approximation has 
been carried out in light cone gauge, A + = 0. This im- 
plies that the + component of the polarization vector e(k) 
is also zero. The physical polarizations of the emitted 
gluon must be transverse to its momentum, e(k) ■ k = 0. 
With these two constraints the two physical polarization 
vectors, i = 1,2, are given by 



-M 



(k) = 0, 



2e ( , 8) • kj 



(9) 



where = (1,0) and = (0, 1) are possible choices. 
For brevity we will suppress the polarization index (i) 
in the following, keeping in mind that the final po- 
larization sum will simply amount to the replacement 

The evaluation of the Feynman diagrams in the GB ap- 
proximation is detailed in |appendix A[ we list the results 
here: 



iM\ 



ii -><?<? g 



2s 



2k x -ei 



iMf ^ qq 9 



l 9 Kj^mn^jl^wsSw's'— (1 — x ) 

q ± k ± 

b 2s 
— i g A^A^„Aj j6 ws 6 w ' s '—2 (1 — x) 

2k ± -ex 



k\ 



iMf ^ qq 9 ~ 
iMV'^ qq ' 9 = 



iMf ^ qq 9 



2s 



9 f a C XilKnn&ws&w' s' — (1 — x) 

9x 

2(qx-kx)-ex 



(«U - kx) 2 



(10) 



Note that the diagrams corresponding to radiation from 
the lower quark line are either exactly (diagram 4) or 
approximately (diagram 3) zero in the A + = gauge. 
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The common factors are already indicative of what 

after squaring the amplitude will be the factorized con- 
tribution of the elastic qq —> qq' process in small angle 
approximation. 

Using [A a , A b ] ;; = if abc ^ii the sum of the five diagrams, 



Eq. [(10)] %M q 

puted 

%M aa' —>aa' a 



i-M-f ^ qq 9 , can be com- 



b 

9 f 1 C ^mnKldws$w's' ~2 (1 



x 2ei 



(11) 



Squaring the amplitude and taking the sums and aver- 
ages over spin states (w, s, w' , s') as well as over the color 
indices yields 



that has also been implemented in the transport model 
BAMPS [ lj. When computing cross sections or rates 
from Eqs. (12) or |(13)| however, the phase space inte- 
gration covers the entire accessible rapidity region and 
factors (1 — x) must not be omitted. Thus Eqs. (12) 
or |(13)| rather than Eq. |(14)| need to be employed when 
computing rates for 2 — > 3 processes in the GB approxi- 
mation. 

The second modification that has been mentioned at 
the beginning of this section is somewhat more subtle. 
In order to arrive at the simple results for the single di- 
agrams as noted in |(10)| when not setting x ~ k_\_/yfs ~ 
0, approximations are necessary that require (see ap- 
pendix A for more details) 



x 2 s > fc 2 



\M 



qq ->qq g\ — Q a „4 \ •*'/ 



+ 



(qU - k±) 2 



or equivalently 
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k + > k y > 



(15) 



(16) 



(12) 



The results Eqs. [(12)1 and EH) 

are thus only valid at 
mid-rapidity, where terms ^x can be neglected, and for 

ing from "the first TcT"thelecondHne"o7[pf is strictly forward emitted S luons ' where x ' s » k l holds ' but J^ 



Note that the algebraic simplification employed when go 



only valid when not screening the propagator terms 1 jk\ 



and l/(q_i_ — kj_) 2 with a Debye mass, cf. Iscction VI 
When introducing the squared amplitude for the elas- 
tic qq' — > qq 1 process in the usual small angle (sa) ap- 
proximation, | Ai qq <^ qq i | 
be rewritten as 



§<? 4 4s 2 M 



Eq. (12 



\M 



qq'->qq'g \ 



12g 2 \M 



qq -*qq 



{l-xf 



ni 



kl(q±-k ± r 



(13) 



explicitly demonstrating that in the high energy limit the 
matrix element of the 2 — > 3 process factorizes into an 
elastic 2 — > 2 part and a gluon emission amplitude. 

The factor (1— x) 2 in Eqs. (12)| or [(l~3"T| lcads to a sizable 
suppression of the amplitude at forward rapidity, w here 
x > k±/y/s, which is immediately evident from Eq. |(4)| 
At the maximal rapidity that is kinematically allowed 
for a given transverse momentum, y max = ]n{k±/'\/s), x 
is 1 and the factor (1 — x) 2 completely suppresses the 
matrix element. While GUNION and Bertsch have the 
(1 — x) terms in their intermediate results for the single 
diagrams, cf. |(10)| they were mostly interested in emis- 
sion spectra in the mid-rapidity region, x ~ k±/^/s <C 1, 
and thus omitted the (1— x) terms from their final results. 
The amplitude then simplifies to 



in the backward rapidity region. The requirement |(15)| 
goes beyond the approximations | (5) | and | (6) | and has not 
been mentioned in the original GB publication. 

Qualitatively, the need for a restriction such as 
Eqs. [(15)1 or p6) can be understood when noting that 
the process qq' — > qq' g is symmetric in interchanging q 
and q' and therefore the resulting amplitude should be 
symmetric in the rapidity of the emitted gluon. Terms 
including x, however, are evidently not symmetric in y, 
cf. Eq. |(4)| and thus the result s [(T2j| and |(13T| do not obey 
the symmetry that should be present. This consideration 
nicely matches the analytic finding that |(12)] and |(13)1 are 
valid only in the forward and mid-rapidity region. 

One can, however, carry out the same calculation in 
the A~ = gauge. In this gauge diagrams 3, 4 and 5 
have sizable contributions, while diagrams 1 and 2 do 
not contribute to the final amplitude. Setting 



(17) 



the final result for the averaged squared matrix element 
in A~ =0 gauge reads 



l-Mqq'^qq'g 



\M 



qq'^qq'g \ 



12g 2 \M 



qq'^qq' \ 



12 9 * \Mqq^qq,\l a {l-x'f 

x 



*i(<n 



*l(qj 



(18) 



(14) 



which is the formula that has been mostly used in the lit- 
erature as the GB matrix element since then and 



This is simply the result obtained from the A + = gauge, 
cf. |(13)| with x being replaced by x' . Since both results 
are also valid at mid-rapidity it is self-evident to combine 
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Egs. [(12)1 and [(18)1 to 

M aa' —¥aa' a 



\2g 2 \M qq >^ qq > ]si 



fci(q ± -k ± )2 



x [(1 - x) 2 e(y) + (1 - x') 2 Q{-y)] 



12g 2 \M qq ,^ qq ,\ sa (l-xy 



fci(q ± -k x ) 2 



where we have defined 



k± 



(19) 



(20) 



In the following we will refer to this result as the im- 
proved GB matrix clement since it is not only valid at 
mid-rapidity, but also at forward and backward rapidity. 
In Sec. [V] we will compare this result to the exact ma- 
trix element without any approximation, which is briefly 
discussed in Sec. IIIII 

The observed factorization of the GB matrix element 
into an elastic part and a gluon emission amplitude in 



the high energy limit, cf. Eq |(13)| immediately implies 
that the GB calculation is also valid for other 2 — > 3 pro- 
cesses, such as qg —> qgg or gg —> ggg. In the high energy 
limit of the GB approximations the specific nature of the 
scattering particles becomes irrelevant for the structure 
of the resulting matrix elements. Therefore, Eq. |(19)| also 
holds for processes other than the one explicitly consid- 
ered here, if one takes into account the different color 
prefactors of the corresponding 2 — > 2 small angle ampli- 
tudes, 



\M 



aa^aa \ 



aa^aa \ 



9V 

! ■ M 



qq'-tqq' 



(21) 



III. EXACT MATRIX ELEMENT 



The exact matrix elements for the 2 — > 3 processes 
involving light partons have been calculated in Refs. 
[24|. For qq' — > qq'g the result is relatively simple, [23 1 



M qq >^ qq 'g\ 2 = y [(S 2 + S' 2 + U 2 + U' 2 ) /it'] [( Pl k) ( P2 k) ( P3 k) (Pik)}- 1 

x {Ci [(u + u') (ss + tt' - uu 1 ) + u(st + s't') + u'(st' + s't)} 

-C 2 [(s + s') (ss' - tt' - uu') + 2tt'(u + u) + 2uu'(t + t')}} , (22) 



with the constants d = (N 2 - I) 2 /AN 3 and C 2 = 
(N 2 - 1)/47V 3 where N = 3 is the number of colors. The 
Mandelstam variables for the 2 — >• 3 process are defined 
as 



s = (pi + p 2 ) 2 t=(pi- p 3 ) 2 u = (pi- pa) 2 
s' = ( P 3+P4) 2 t' = (p 2 - Pi ) 2 u' = (p 2 -p 3 ) 2 (23) 



However, the expression for the matrix element has been 
simplified to such an extent that it is not obvious any- 
more how to identify the internal propagators which are 
usually screened with a mass of the order of the Dcbye 
mass in thermal QCD. We do not explicitly quote the 
result for qg — > qgg here, it is also given in [23[. 

The expression for gg — > ggg is considerably more com- 
plicated, but can still be ex pre ssed in a relatively compact 
form due to its symmetry |23| 



,2 " 6 



\M ag ^ gag | = — [N 3 /(N 2 - 1)] [(12345) + (12354) + (12435) + (12453) + (12534) 

+ (12543) + (13245) + (13254) + (13425) + (13524) + (14235) + (14325)] 
[(P1P2) 4 + (P1P3) 4 + (P1P4) 4 + (P1P5) 4 + (P2P3) 4 ] 

(PlP2 ) (PlP3 ) (PlPi) (PlPh) (P2P3 ) (P2P4 ) (P2P5 ) (P3P4 ) {P3P5 ) (P4P5 ) 

_! [(P2P4) 4 + (P2P5) 4 + (P3P4) 4 + (P3P5) 4 + (paPb) 4 } ^ 

(PlP2 ) (PlP3 ) ipiPi )(PlPb) (P2P3 ) {P2PA ) (P2Ph ) (P3P4) (P3P5 ) (P4P5 ) 



with (ijklm) = (piPj)(pjPk)(PkPi)(piPm)(PmPi) and p§ notation p§ = k. As for the other matrix elements, it 
being the momentum of the third outgoing gluon. In our 
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is not a priori obvious how one would screen this matrix 
element. For the purpose of this work we introduce a cut- 
off prescription in Sec. [V] that regulates the divergencies 
when calculating the total cross section and allows for 
consistent comparisons to the GB result. 



IV. PHASE SPACE INTEGRATION FOR CROSS 
SECTIONS 

In the previous sections the leading order pQCD ma- 
trix elements for radiative processes have been introduced 
both in their exact form and in the (improved) approx- 
imation by Gunion and Bertsch. These matrix ele- 
ments depend on four independent phase space variables. 
A possible choice for these variables that is most appro- 
priate for the GB approximation of the matrix element is 
shown later in this section, cf. Eq. |(26)| In the following 
we outline the procedure for the phase space integration 
that is needed to obtain the total cross section from the 
matrix elements presented in Sections [TT| and HID Full 
energy and momentum conservation are taken into ac- 
count. In general, differential cross sections (e.g. da/dy) 
can be obtained from this procedure by projecting out 
variables with delta functions (e.g. 8(y)). It is these dif- 
ferential and total cross sections that are then used in 
Iscction Vl to numerically check the improved GB approx- 
imation against the known exact result as they are easier 
to visualize than the bare matrix elements. 

The total cross section of a 2 — > 3 process is defined by 



02->3 



d 3 p 3 



d 3 p 4 



d 3 fc 



2£i 2E 2 Wrci v J {2n) 3 2E 3 (27r) 3 2^ 4 (27r) 3 2fc° 



x (27r) 4 d {4) ( Pl +p 2 - (p 3 +P4 + k)) \M 2 ^ 3 \ 



(25) 



where we used the same notati on as before. The r elative 
velocity is defined by u rc i = \J (-Pf i^) 2 — mfm 2 / 'E\E 2 
and the symmetry factor v — n\ accounts for n identical 
particles in the final state. Note that this symmetry fac- 
tor is not needed when computing cross sections from the 
(improved) GB matrix element. In this case the specific 
choice and identification of the outgoing momenta, where 
P5 = k is the momentum of the radiated gluon, selects 
a specific configuration and thus obviates the need for a 
symmetry factor 1/v. This is taken into account for all 
calculations involving the original or improved GB ma- 
trix elements in this paper and has also been correctly 
implemented in BAMPS from the beginning [lij . 

Equation |(25)| can be used to calculate the total cross 
section of the GB and exact matrix elements by perform- 
ing the nine-dimensional integration. However, it is nu- 
merically more efficient to analytically integrate out the 
delta function to end up with only a five dimensional inte- 
gration, which can be reduced to a four dimensional inte- 
gration by performing one angle integration directly. Af- 
ter this and switching to the GB coordinates from Sec.HTl 



the total cross section reads \\A 



S/4 g/4 l/max 

a ^ 3 = 2Ehl- s J dql S dkl J dy 

J„in 



x \M 



2^3 



n — 



F=0 



(26) 



where 7/3 denotes the rapidity of particle 3 and <j) the angle 
between qj_ and kj_ . The available phase space limits the 
rapidity to y max /min = i arccosh \yfs/2k±]. Due to the 
transformation of the delta function, (dF/dy 3 ) 1 must 
be summed over the roots of 

F = (pi + p 2 - p 3 - k) 2 

= s — 2^/s (q± cosh 1/3 + k± coshy) + 2q±k± cos<f> 
+ 2q± k± (cosh y 3 cosh y — sinh y 3 sinh y) . (27) 



We have numerically checked that both Eqs. | (25)] and 
|(26)| give the same results for the GB and exact matrix 
elements, respectively. 



V. NUMERICAL COMPARISON 

We have carried out extensive numerical comparisons 
between the exact matrix element and the matrix ele- 
ment calculated in the GB approximation. In order to 
obtain a picture as complete as possible we have inves- 
tigated differential cross sections such as da/dy, da/dx, 
da/dk±dq± in addition to the total cross section. All 
calculations take into account the full kinematics accord- 
ing to Eqs. I (2571 or I (26)] ensuring energy and momentum 
conservation. 

Both the exact and GB matrix elements are diver- 
gent for infrared and collinear configurations. In thermal 
quantum field theory these divergencies can be cured by 
loop resummations of the propagator, which leads to an 
extra term with the self energy in the propagator. Phe- 
nomenologically, this self energy can be mimicked with a 
Debye mass m d that modifies the propagator terms and 
makes them infrared-safe, 



1 



1 



(28) 



D 



However, due to their compact notation it is not 
straightforward to identify the propagator terms in the 
exact matrix elements. Therefore, for all numerical cal- 
culations of total or differential cross sections from both 
the exact and the GB matrix elements we choose a sim- 
ple cut-off procedure to cure the divergencies. The in- 
tegrand is set to zero when the scalar product of any 
incoming or outgoing 4-momenta is smaller than a cut- 
off A 2 , pi ■ pj < A 2 with i, j = 1...5. Formally this is 
expressed by multiplying the integrand in Eq. |(25)| with 
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Figure 2. Differential cross section da/dy for the process 
qq' — > qq' g calcu lated with the exact (Eq. (22) I, GB (14) 
GB with (1 - 



(12) and improved GB with (1 — a;) 
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matrix element. 

Q(pi ■ Pj — A 2 ), where A 2 is chosen to be proportional to 
the Dcbye mass, 



A 2 = 



(29) 



It is immediately obvious from Eq. |(24)| that this prevents 
divergencies of the integrand. For propagator terms this 
scheme acts similar to the Debye screening, but the re- 
strictions to the phase space potentially go beyond Dcbye 
screening only propagator terms. The resulting numeri- 
cal values for the total or differential cross sections need 
thus not be the physically correct values and are conse- 
quently given in arbitrary units where necessary. Still 
this cut-off prescription allows for consistent and well- 
defined comparisons between the GB approximation and 
the exact leading order matrix elements which is the fo- 
cus of this study. If not mentioned otherwise, we set 
e = 0.001 to reduce the screening effect as much as possi- 
ble. In the last part of the section we will discuss the im- 
pact of different e and compare this cut-off procedure to 
the standard Debye screening procedure for more physi- 
cal scenarios. 

The calculations in the remainder of the section are 
done for a temperature of T = 400 McV. The coupling 
is set constant, a s = 0.3. If not stated otherwise, we use 
the average thermal value for massless particles for the 
squared center of mass energy, s = 18T 2 ~ 2.88 GeV 2 , 
and determine A from the usual gluon Debye mass for 
Boltzmann statistics, m 2 D = ^-{N c + nf)T 2 . For nj = 3 
at the given temperature and coupling, the Debye mass 



m 2 D ~ 0.73 GeV 2 . 



Figurc[5]compares the rapidity spectrum of the emitted 
gluon as given by different approximations of the matrix 
element by depicting the differential cross section da/dy 
for the process qq' — > qq' g. This process is the clearest to 
study since only one (the emitted) gluon is involved. Fur- 
thermore, it is the process that has been studied in the 
original GB publication, cf. Isection III The plot nicely 



demonstrates the shortcomings of the standard GB ma- 
trix element, Eq. |(14)[ in the forward and backward re- 
gion. 

The exact matrix element, Eq. |(22)| only peaks at mid- 
rapidity and does not have any sizable contribution at 
forward or backward rapidity. The gluon emission into 
the forward and backward region is suppressed since the 
phase space is occupied by the two quarks which just 
scatter off with a small angle. Furthermore, the cross 
section is symmetric in the rapidity of the emitted gluon 
as it must be for this process. The standard GB ma- 
trix element without the (1 — x) 2 term (Eq. |(14)[ ) is very 
similar at mid-rapidity, but has two additional large con- 
tributions at forward and backward rapidity. Thus this 
matrix element would allow that, for instance, the gluon 
is emitted into the backward region and that one of the 
quarks is located at mid-rapidity. Although the matrix 
element itself is symmetric in y, the curve is not since in 
our implementation the p z of particle 3 is restricted to 
positive values to ensure that the particle only scatters 
off with a small angle Q 

The matrix element with the GB approximation and 
the term (1 — x) 2 , Eq. |(12)| has the same value at mid- 
rapidity as the other matrix elements since x is small in 
this region. However, at forward rapidity the (1 — a;) 2 
factor leads to a significant reduction compared to the 
pure GB result. Here the curve lies right on top of the 
curve from the exact matrix clement. Nevertheless, as 
discussed in Sec. [ill the curve with (1—x) 2 is not symmet- 
ric anymore due to the asymmetric y-dependence con- 
tained in x. The last curve in Fig. [5] shows the cross 
section of our proposed improvement of the GB matrix 
element, Eq. |(19)| At forward and mid-rapidity it agrees 
with the previously mentioned curve. At backward ra- 
pidity the modified (1 — x) 2 term removes the excess and 
reconciles the shape of the GB approximation with the 
symmetry of the process. The overall agreement between 
the differential cross sections from the improved matrix 
element and from the exact result is remarkably good. 
The difference between the approximation and the ex- 
act result is on the level of a few percent over the entire 
y-rangc. This striking agreement is also visible in other 
differential cross sections such as da/dq^dk 2 ^ and da/dx 
(not shown here). 

Figure[5]was done for the simplest process, qq' — > qq'g, 
where only one gluon is involved. However, the findings 
also hold for qg — > qgg and gg — !> ggg. Since gluons are 
indistinguishable particles we plot the da/dy for qg — > 
qgg of all outgoing gluons in Figure [31 Since we choose 
the incoming gluon to be particle 1 and the process favors 



1 This constraint is necessary if one transforms from Mandelstam 
t to for such a small angle favoring process, since q 2 ^ can be 
small for a outgoing particle which moves exactly in the other 
direction than it moved before the reaction, but in this case the 
scattering angle is not small anymore, thus, t becomes large and 
t ~ q 2 ^ does not hold anymore. 
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Figure 3. Differential cross section da/dy for all outgoing 
gluons for the process qg — s> qgg calculated with the exact 
and improved GB matrix element. 
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Figure 4. Ratio of total cross section <72->3 of the improved 
GB matrix element (Eq. | ( 19) [ ) to the exact one (see Sec. 
as a function of s for qq — > qq'g, qg — > qgg, and gg —>■ ggg. 
In addition, the ratio of the original GB cross section (14) as 



it was previously implemented in BAMPS to the exact cross 
section is plotted for gg — > ggg. 



small angle scattering, there is a peak at forward rapidity 
in addition to the gluon peak at mid-rapidity. Again the 
curve from the improved GB matrix element reproduces 
the curve from the exact one very well. This also holds 
for gg — > ggg, where one gets an additional gluon peak 
in the backward region which corresponds to the third 
outgoing gluon. 

Since all differential cross sections for the improved GB 
and exact matrix elements agree nicely, it is not a sur- 
prise that this is also true for the total cross sections. In 
Figure 2] the ratios of the total cross sections <72->3 for 
different processes from the improved GB to those ob- 
tained from the exact matrix elements are depicted as 
a function of the squared center of mass energy s. The 
exact and improved cross sections agree very well for all 



Figure 5. Ratio of the total cross section a of the improved GB 
matrix element to the exact result for qq' — s> qq'g, qg — > qgg, 
and gg — > ggg as a function of the screening cut-off e from 
Eq. (29) In addition, the ratio of the original GB cross section 
of gg — > ggg to the exact is depicted as it was previously 
implemented in BAMPS. 



processes and virtually all s. Only at very small s, below 
the thermal average s = 18T 2 ~ 2.88 GeV 2 , there is a 
slight discrepancy up to about 40%. The original GB 
cross section for gg — > ggg as it was previously imple- 
mented in BAMPS is about a factor of three larger than 
the exact cross section. The same is true for the other 
processes (not plotted). 

Figure [5] shows the dependence of the total cross sec- 
tion on the cut-off e = A 2 /m 2 D . The improved GB matrix 
element works best for small e where the matrix elements 
arc dominated by small k± and q_\_. At larger e most of 
the phase space at small k± and q± is cut away and other 
phase space regions play the dominant role, in which the 
GB approximation is not as good. Hence, the ratio rises. 
At e > 0.6 (indicated by the grey band in the figure) the 
cut-off is so large that the entire available phase space 
is cut away and the GB as well as exact cross sections 
are zero. For such severe cut-offs the comparison scheme 
is not reliable any more and the previously introduced 
Debye screening would be more realistic. However, the 
comparison to the exact matrix element with standard 
Debye screening is not done in this paper since the iden- 
tification of the propagators in the compact results of the 
exact matrix elements is not known. Also note that Fig. [5] 
is made for thermal s = 18T 2 . For larger values of s the 
cut-off parameter e at which the ratio deviates from unity 
becomes significantly larger. Thus, in summary, the im- 
proved GB approximation works best for larger s and/or 
smaller screening cut-offs e. 

Our improved GB result is therefore an especially good 
approximation for studying jet effects in heavy-ion colli- 
sions. For bulk dynamics, where the mean s is smaller, 
it is important how strong the screening effects actually 
are and if one is already in the regime of large e where 
the improved GB matrix clement might deviate from the 
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Figure 6. Differential cross section da/dy for the process 
qq' — > qq'g calculated with the improved GB matrix element 
with the standard Debye screening and the cut-off procedure 
for different e. 



Figure 7. As Fig. [5] but in addition the improved GB cross 
section with the exact 2 — > 2 matrix element (22 t) instead 
of the 2 — > 2 matrix element in small angle approximation 
(22 q ± ). 



exact matrix element. To explore this in more detail, in 
Fig. [6] we explicitly compare the improved GB differen- 
tial cross section obtained from the cut-off procedure to 
the one obtained using the standard Debye screening ac- 
cording to Eq. |(28)| which is for qq' — > qq'g (cf. Eqs. (12) 
andpjjflg 



\M, 



i 2 o) 2 



kl 



D 



(q± - kj.)s 



D 



(30) 



Now to get a feel for the quality of the improved GB 
approximation for Debye screened thermal processes one 
can either try to match the total cross sections or the 
differential cross sections da/dy at mid-rapidity. Doing 
this, the necessary parameter e can be determined to be 
of the order of e « 0.4 to 0.5. As can be seen from Fig. [5] 
the deviations between the improved GB and exact result 
in this e region vary roughly between cr/cr xact ~ 4 and 
c/cexact ~ 9. However, Fig. [5] also clearly demonstrates 
the distortion of the phase space due to extreme choices 
of the cut-off. Large rapidities arc cut away entirely 
while the standard Debye screening procedure leaves the 
whole phase space available. It is questionable whether 
the strong cuts in the cut-off procedure make physical 
sense. The comparison is thus a rather qualitative one 
and should be regarded as a first estimate of the qual- 
ity of the (improved) GB approximation in the region of 
thermal processes with a standard Debye screening. 



2 The implementation in BAMPS does not screen the term in 

the bracket since small k± are regularised by the LPM effect [25l . 
However, since the LPM effect is not employed in this study, we 
regularise also this term with the Debye mass. 



As discussed in Isection III the factorization of the GB 
approximation into an elastic 2 — > 2 part and a radiation 
amplitude features the elastic amplitude in small angle 
approximation, |A^2-s.2| sa oc l/q\. If one, however, in- 
stead employs the exact binary matrix element for the t 
channel which is 



Mo 



I 2 

>2 OC 



1 



(31) 



the agreement between the exact and the improved GB 
result becomes much better at large e, which is illustrated 
in Figure [7] The ratio is of the order of 1 over the entire e 
range and the deviation is at most 50 % for all processes. 
In contrast to the previous result with the approximated 
2—^2 matrix element (cf. Fig. [5]), the ratio for small e is 
not 1 anymore, but 0.9. This is probably due to the fact 
that different approximations are employed for the 2 — > 2 
and the radiative part of the matrix element. However, 
since the numerical agreement is so good also in the re- 
gion of thermal processes at larger e, the combination 
of the improved GB approximation |(19)| for the radia- 
tion amplitude with the exact 2—^2 part |(31)| might be 
implemented in our transport model BAMPS. 



VI. SUMMARY 

The inelastic gluon Brcmsstrahlung matrix element de- 
rived by GUNION and Bertsch for qq' — s- qq'g has been 
compared to the exact matrix element for this process 
by performing full phase space integrations and analyz- 
ing the total and differential cross sections. Although it 
agrees very well with the exact result at mid-rapidity (the 
region in which Gunion and Bertsch were mainly inter- 
ested in) large deviations have been found at forward and 
backward rapidity. Based on a detailed analytic investi- 
gation of the underlying approximations, we propose an 
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improved version of the GB matrix element which agrees 
very well in all phase space regions with the exact result. 
This remarkable agreement holds for virtually all center 
of mass energies y/s and most cut-off parameters e as 
well as for the other radiative processes qg —> qgg and 
99 -> 999- 

In contrast to the claim in Ref. [2l[, the deviation 
between the GB matrix element and the exact result 
does not originate from a miscounting of symmetry fac- 
tors, but lies deep in the approximations of GB in for- 
ward and backward rapidities. Although the reasoning 
in Ref. [2l[ accidentally (and thus only effectively) holds 
for gg — > ggg due to the symmetry of the process, it fails 
for asymmetric processes such as qg — > qgg or qq' — > qq'g. 

The implementation of the new improved GB matrix 
element in our transport model BAMPS is currently un- 
der way. From Fig. 2] we expect at most a factor of 3 
difference in the rates, although this should be reduced 
by the iterative computation of cross sections due to the 
implementation of the LPM effect [25| . The difference in 
the phase space sampling of the outgoing particles will 
have a large impact in particular on jet obscrvables. For 
instance, we expect much less suppression of high energy 
particles, which will increase the small nuclear modifi- 
cation factor of BAMPS [13, HE| and bring it closer to 
the data. The effect on bulk observables will also be 
of great interest. Although we expect a decrease of the 
ellip tic flow and an increase of the shear viscosity 
[20|, one must wait for the numerical results to see how 
large the effect will actually be. As part of this study we 
will also implement the running coupling for elastic and 
radiative processes which is expected to lead to a rising 
Raa curve with transverse momentum. Furthermore, we 
plan to generalize the new improved GB matrix element 
also to the heavy flavor sector and study its impact on 
the heavy quark energy loss [26l - [28j j . 
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Appendix A: Detailed GB calculation 

The purpose of this appendix is to provide some more 
details on the computation of the qq' — > qq'g matrix el- 
ement in the Gunion and Bertsch approximation, cf. 
Sec. [H] and the diagrams shown in | Figure l) The focus 
is on highlighting crucial steps and approximations that 
are needed to arrive at the results in Eq. |(10)| 

We define 



(Al) 



where s and r are the spin indices of the spinors u and u. 
Setting p' = p + q and assuming that q <C p, the Gordon 
identity [H] for pi)] can be simplified to give 



J? s (p,p + q)~(2p + qr6 r 



(A2) 



This eikonal approximation essentially amounts to de- 
scribing the problem in scalar QCD and is responsible 
for the fact that the amplitudes of different 2—^3 pro- 
cesses only differ by the color factor of the factorized 
elastic amplitude, cf. Eq. |(21)| Furthermore we define 

T%(p, k,p') = u r (p) r (* + m) 7 " u s (p') 

= £j£(p,*)J£(fc,j/) • (A3) 
t 

With these definitions and the kinematics as given in 
Sec. |TT] we can compute the matrix elements of the dia- 
grams. 



Diagram 1 



The Fcynman rules for diagram 1, cf. |Figure 1[ give 

<'a Fl 



q% , + k y rSs(P3,P3 + k,Pl) 

■F\ S WS S W ' s ' , \u fn \ 

— — 2 (2pi + q )" {2p 2 - g) A 

x (2 Pl +2q-kre* a (k) 



(A4) 



with F\ = —i g 3 A°,- A^ w A^ ac . Employing the kinematic 
approximations from Sec. [H] to the single terms 

2 (k ± - xq_L) 2 



(Pi + qY 

{2p l + qY{2p 2 -q) ll 
(2 Pl + 2q- k) a e* a 



(1 — x)x 



2s 



x 



(A5) 
(A6) 
(A7) 



together with q ~ — qj_ and |(6)| yields the result listed 



(15) 


in 


and 


(6) 



in addition to the standard GB approximations (5) 



10 



Diagram 2 
The Feynman rules for diagram 2 give 



K a s (P3,Pi - k, Pl ) 



q 2 {p 1 - k) 
F 2 5 WS 8. 



q 2 {pi - k) 2 
x (2pi - fc)X(fc) 



(2p 1 + g-2*)"(2 J j2-g) M 
(A8) 



with F 2 = —ig i \%Xj. nn \ b ii 5 ac . The kinematic approxi- 



mations 



A- 2 



(A9) 



(2pi + g - 2fc) A '(2p 2 - g)„ ~ 2s (1 - x) (A10) 



(2 Pl -fc) Q e* ~-k ± -ej 



(All) 



together with g 2 ~ — g 2 le ad to the result listed in 
Eq. |(10)| Again Eq. |(A10)| requires the restriction to 
forward emission, x 2 s 3> k ± , in addition to the standard 
GB approximations |(5)| and |(6)| 



Diagram 3 

In order to explicitly compute diagrams 3 and 4 the 
components of the momentum transfer q need to be rede- 
termined from the on-shell conditions (p 2 = (pi + q) 2 = 
and p\ = (p 2 — q — k) 2 = 0) and the choice of keeping k 
in the form given in Eq. |(7)| 



-Xy/S, -=,qj 



(A12) 



With this definition diagram 3 can be evaluated as 
F 3 



,P4 + k,p 2 )e* a (i 
(2p 1 +qY{2p 2 -q) ll 



q 2 (pi + k) 
xT a J s ,(p^p 4 + k,p 2 )e* a (k) 

$ws$w' s 



q 2 (P2 - q) 2 

x(2p 2 -2q-k) a e* a (k) (A13) 

with F$ = —ig 3 X^i X c mj ,X h j n 8 ac . Employing the usual kine- 
matic approximations simplifies the single terms 

(p 2 - q) 2 ~ xs (A14) 

{2p l +qY{2p 2 -q)^s{2~x) (A15) 

(2p 2 -2 g -A : ) Q e;~2(q ± +k ± )- e± . (A16) 



Combining these approximations yields 

iM f^B ~ F 3 <WW^ — (q±+k ± )-e ± (A17) 

Qx s 

which is suppressed by 1/s compared to the other dia- 
grams listed in |(10)| and thus in the GB approximations 

justifies setting iM. qq ~* qq 9 ~ 0. 



Diagram 4 



Diagram 4 is given by 



iMf ^ qq 9 



q 2 (p 2 - k) 2 
x (2p 2 -k) a e* a (k) 



J% s (P3,Pl)9lxv 

,P2 ~ k,p 2 )e* a (i 
{2p 1 +qY{2p 2 -q-2k) ti 



q 2 (p 2 - k) 2 
x TZ> s >{paiP2 - k,p 2 )e* a {k) 

^4 $ws$w' s' 



(Ais) 



with F 4 = -ig 3 Xl\ c jn X h mj 5 ac . Since (2p 2 - k) a e* a = 
without any further approximation, diagram 4 does not 
contribute. 



Diagram 5 

Diagram 5 that describes the gluon emission off the 
exchanged gluon is the most tedious one to evaluate. It 
is given by 



iMf ^ qq 9 



Jws(P3,Pl)Jw' S '(P4,P2) 



q 2 {q — k) 2 
xe* a (k) [g^{k-2qr + g va {q + kY 
+g a »(q-2k) v ] 

^^W [A + B + C] (A19) 

with F 5 = g 3 f abc X a a X c mn and 

A = [(2pi + q - kY(q + k)„} [{2p 2 - qYC] 

cl—^^q ± -e ± (A20) 
1 — x x 

B = [{2p 2 - qf{q - 2%] [(2 Pl + q - k) a e* a ] 
~-4sk_L-e_L (A21) 

C = {(2 Pl + q - kY(2p 2 - g)„] [{k - 2q) a e* a ] 
~ 2s(2-x)q_L-e_L . (A22) 

The term A is not proportional to s and thus can be 
neglected in respect to B and C. Together with the ap- 



proximations for the propagators, q 2 ~ — q 2 and 



(q-kY 



(A23) 



(l-x) ' 

this leads to t he r esult in Eq. |(10)| where additionally 
approximation [(6)j needs to be employed in collecting the 
above given terms. The constraint x 2 s 3> k 2 ^ is needed in 
the computation of diagram 5 when approximating the 
first term of C as given in |(A22)| 
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